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We present a mapping which associates pure JV-qubit states with a polynomial. The roots of the 
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I. INTRODUCTION 

Considerable effort is spent in developing methods for 
the detection and classification of entangled states. One 
important aim is to find ways to detect the separabil- 
ity of mixed states consisting of an arbitrary number of 
subsystems. While a general, easily computable, method 
to detect the separability of arbitrary mixed multipartite 
states is still lacking, some partial results exist. Maybe 
the most famous separability condition for mixed states 
is the positive partial transposition, also known as Peres- 
Horodecki criterion This method is simple and easy 

to apply, but it can be used to detect only bipartite sepa- 
rability. Therefore various separability conditions which 
work in an JV-party setting have been developed. Exam- 
ples of these are permutation criteria, where the indices of 
the density matrix are permuted [3j , the use of quadratic 
Bell- type inequalities [3] , algorithmic approaches @ , and 
the use of positive maps [6j. For a more comprehensive 
list, see @, H|. In the case of pure states the situation 
is simpler. A pure A-partite state is separable if and 
only if all the reduced density matrices of the elemen- 
tary subsystems describe pure states. Alternatively, in a 
bipartite case, separability can be determined by calcu- 
lating the Schmidt decomposition of the state. Unfortu- 
nately, the concept of the Schmidt decomposition cannot 
be straightforwardly generalized to the case of A separate 
subsystems d, [l(| ■ In addition to these two well-known 
methods, various other approaches to the pure state sep- 
arability have been discussed. A separability condition 
based on comparing the amplitudes and phases of the 
components of the state has been discussed in [IJ [l2j ■ 
It has been shown that the separability of pure three- 
qubit states can be detected by studying two-qubit den- 
sity operators [l3| and expectation values of spin opera- 
tors 0, . Separability tests based on studying matri- 
ces constructed from the components of the state vector, 
known as coefficient matrices, have gained attention re- 
cently [H-[l|. 

In this article we present a mapping which associates 
the pure states of an A-qubit system with a polynomial. 
The roots of the polynomial determine the state com- 
pletely and vice versa. We show that this polynomial 
establishes a simple way to test the separability of pure 
JV-qubit states and to study the number of unentangled 
particles. The idea to associate a state of a quantum 



mechanical system with a polynomial is not new. Al- 
ready in 1932 E. Majorana presented a polynomial, nowa- 
days known as the Majorana polynomial, which he used 
to show that the states of a spin-5 particle can be ex- 
pressed as a su perp osition of symmetrized states of 25 
spin-i systems [l9l |20|. This decomposition, the Majo- 
rana representation, has been relatively unknown for a 
long time. However, it has recently found applications 
in may different fields, such as in studying the symme- 
tries of spinor Bose-Einstein condensates |2ll - l24| , in the 
context of reference frame ali gnm ent (25j |. in helping to 
define anticoherent spin states [261 ] , and in calculating the 
spectrum of the Lipkin-Meshkov-Glick model [13,[28|. It 
has also been used to give a graphical representation for 
the states of an n- level system |29j . 



The states of an A-qubit quantum register can be 
viewed as the spin states of a particle with spin 5 = 
(2 N — l)/2. Therefore, expressing the pure states of an 
A-qubit system utilizing the Majorana representation re- 
quires the use of 2 N — 1 spin-i systems. In the approach 
we present in this article only A two-level systems are 
needed to characterize the states of this system. The 
Majorana representation is useful in studying the behav- 
ior of spin states under spin rotations as a spin rotation 
of a spin-5 particle is equivalent with rotating the states 
of the constituent spin-^ particles However, when 

discussing the states of an A-qubit quantum register, this 
property is not very helpful and therefore the benefits of 
the Majorana representation cannot fully be taken advan- 
tage of. In this case the simplified description presented 
in this article becomes useful. 



This article is organized as follows. In Sec. II we intro- 
duce a mapping between the pure states of an A-qubit 
quantum register and polynomials. We argue that the 
roots of a polynomial determine a unique state and vice 
versa. In Sec. Ill we calculate the polynomial of sep- 
arable pure states and derive a necessary and sufficient 
condition for the separability of an arbitrary pure A- 
qubit state. We also briefly discuss the generalization of 
the polynomial approach to systems containing JV copies 
of an ft,-level system. In Sec. IV we show how the poly- 
nomial can be used to study the number of unentangled 
qubits. In Sec. V we present the conclusions. 
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II. CHARACTERISTIC POLYNOMIAL 

We denote the basis of the qubit j by {\Q)j, so the 

basis vectors of an iV-qubit quantum register can be cho- 
sen as |io*i • • • iN-i) = Ko)o <8 \h)i ® • • • ® |«iv-i)jv-i, 
where every ij £ {0,1}. Each natural number < 
i < 2 N — 1 can be written using binary notation as 
i = Y^j^o i where ij £ {0,1}. Using this we can 
associate the basis vector |zo*i • ■ -in-i} with \i)d- Here 
the subscript d shows that decimal notation is used to 
label the basis states. Let 

2 N -1 

cb=Y, Ci\i) d , Ci£C, (1) 

i=0 

be some, possibly unnormalized, state vector of an N- 
qubit system. We associate this vector with the polyno- 
mial 

2 N -1 

P{4>;x) = ( 2 ) 
i=0 

which we call the characteristic polynomial of (j>. By the 
fundamental theorem of algebra, this polynomial can be 
written in a unique way as 

fc-i 

P(4r,x) = C k H(x-Xj), (3) 

3=0 

where {xj \j = 0, 1, . . . , k — 1} are the roots and k is the 
degree of P(</>; a;). If k — we define FJ~^ (a: — Xj) = 1. 
The set of vectors {ccf>\c £ C,c ^ 0} determines a unique 
set of roots and each set of roots {xq, x\, . . . , Xk-i} deter- 
mines the vector <p up to normalization and phase. There- 
fore we have a bijective map between the pure states of 
an ./V-qubit quantum register and the roots of complex 
polynomials of degree k < 2 N — 1 [3l|. Explicitly, the 
components of <p are determined by the roots through 
the formula 

Q = ( — 1) y t X jo X jl X 32 ' ' ' X jk-l-i' (4) 

jtt<ji<h<---<jk-i-i 

where i = 0, 1, 2, . . . , k — 1 and we have chosen Cfc = 1. 
The roots contain the same amount of information on 
the system as the state vector (j). In particular, all the 
entanglement properties of (f> are encoded in the set of 
roots corresponding to (j). With the help of the roots 
the state <j) can be given a geometrical representation as 
2 N — 1 points on the Bloch sphere, see Ref. (30j . 



III. SEPARABLE PURE iV-QUBIT STATES 

In this section we show how the separability of <j> can be 
detected with the help of P(<j>;x). In order to do so, we 



first calculate the characteristic polynomial of a separable 
state. Any separable pure state <p s can be written as 



JV-l 




3=0 
N-l 



= ®fo|0)i + &i|l)i) eC. (5) 

3=0 

Assume that \l)d is a basis state of an L-qubit system 
and that \m}d is that of an independent M -qubit system. 
Using the binary expressions for / and m it is easy to see 
that 

\l)d\m) d = \l + 2 L m) d (6) 

holds for the tensor product of \l) d and \m)d- Here and in 
what follows we omit the tensor product symbol. Let £ L 
and £ M be states of L-qubit and Af-qubit quantum regis- 
ters, respectively. Then we can write £ L = Yli=o 1 

and £ M = Y^Lq 1 d' \i')d- K G (C 2 ) L+M can be written 
as <f> = then 

2 L -12 M -1 

( t> = E E ^\iWU (7) 

i=0 i'=0 
2 i -12 M -l 

= E E ^1+^1*, (8) 

1=0 i>=0 

where we have used Eq. ©. Consequently, the charac- 
teristic polynomial of <j> becomes 

2 L ~l 2 A/ -1 

j=0 i'=0 
2 L -1 2 M -1 

= (E^ 1 ) E ^^ L r 

i=0 i'=0 

= p{ev x )p{e i i x2L ). (9) 

Therefore, if the state of the quantum register is the 
product of an L-qubit state and an M-qubit state, the 
characteristic polynomial factorizes as the product of the 
polynomials of the two states. In the polynomial of the 
M-qubit state the variable x is replaced by x 2 . Using 
Eq. ((9]) it is easy to calculate the characteristic poly- 
nomial P(4> s \ x) of a separable state <fi s = 4>o4>i ■ ' 1 4>n-i 
given by Eq. ©. By defining <pj. N = <f>j<j>j+i ■ ■ ■ 4>n-i, 
so that 4>j- t N = 0j0j_)_i ; jv, and using Eq. © repeatedly 
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we get 

P(<f> s ;x) =P(fo,x)P(fa ]N ,x 2 ) 



= P{^ ,x)P^ 1 ,x 2 )P(^ 2 . N ,x 4 ) 

= P{fo,x)P(0 1 ,X 2 )P{fo 1 X 4 )P(fo; N ,X 8 ) 

N-l 

= WP^^x 23 ) 

3=0 

N-l 

1 1 ' ' • 'v 2 >• (io) 

3=0 

We see that the characteristic polynomial of a separable 
state can always be written in the form of (|10[) . On the 
other hand, there always exists a separable state whose 
characteristic polynomial is given by Eq. (|10[) . namely 
the state <j> s . From the definition of P(<fi; x) it follows that 
if P(<f>;x) = P((j>;x), then necessarily <fi = (j). Therefore 
(j) s is the unique vector which gives rise to the polynomial 
of Eq. dTU|) . In conclusion, a pure iV-qubit state (f> is 
separable if and only if P((f>; x) can be written as in Eq. 
(ITU1) . The roots of this equation are 



1/2 3 



0,1 2< !. (11) 



where bj has to be nonzero. If bj is zero the degree of the 
polynomial is decreased by 2 J from the maximal degree 
2^-1. 

The separability of a state <f> can be determined by 
calculating the roots of P((f>, x) and checking if they are 
of the form given by Eq. (fTTj) . These calculations can 
in practice turn out to be very complicated. It may be 
computationally demanding to achieve accurate enough 
results in order to reliably see how the roots are dis- 
tributed in the complex plane. This is partly related 
to the fact the degree of the polynomial P(<fi; x) can be 
2 N — 1, which grows rapidly with N, rendering the calcu- 
lation of roots time-consuming for large N. However, we 
will show next that the roots of P(<fi; x) can be expressed 
in a simple way in terms of the components {Ci} of the 
state vector if is separable. Let 4> s be the separable 
state given by Eq. ([5]). When this vector is written in 

2 N — 1 

the form <f) s = 5Ji=o Cilfydi the components Ci are eas- 
ily obtained by noting that ij = (ij = 1) corresponds 
to aj (bj): 



N-l 

ma 

3=0 



h) a 3 



(12) 



Here we have used the binary form of i, that is, we have 
written i = Xl^Jo 1 ■ We assume that C k ^ 0, C k +i = 
■■■ = C 2 «_ 1 = 0, so that the degree of P((/> s ;x) is k. 



By writing k = X^jLo 1 we see that if kj = 1, then 
(k - 2 ] )i = ki - Sji, I = 0, 1, . . . , N - 1. Using this and 



Eq. (fT2j) it is easy to see that now aj/bj — C fc _ 2 j /Ck- 
On the other hand, if kj = 0, then (k + 2 J )i = ki+Sji and 
Eq. (fT2"j) gives bj/aj = Cfc +2 j/Cfc = 0. Summarizing, 



c 



fc-2J 

bj ~ C k 

h = o 



if kj = 1, 
if kj = 0. 



(13) 



Using Eq. (jTTJ) we immediately see that the k roots of 
P(4> s \x) are 



c k 



e 2J 



m = 0, 1, 



(14) 

where j takes those values for which k } ■ = 1. On the other 
hand, if the roots and their multiplicities are known, the 
polynomial can be determined up to a multiplying con- 
stant. In particular, if x = is a root, then its multiplic- 
ity has to be equal to the lowest power of the polynomial. 
In conclusion, an arbitrary pure state <f> is separable if and 
only if 

(la) All the numbers Xj m given by Eq. (|14p are 

roots of P(4>; x). 
(Ib) The number of Xj m equaling zero is equal to 

the lowest power of P((j);x). 

An alternative formulation is that (f> is separable if and 
only if the quantity 



E E \ P (^ X 3m)\ 



(15) 



j, kj—1 m=0 



equals zero and Condition (Ib) holds. Note that if k = 
the state is separable. 

If a state is found to be separable, then the one-particle 
states it consists of can be explicitly constructed with 
the help of the ratios aj/bj given by Eq. (TIB")) . We now 
present some examples of the detection of separability of 
states with several freely varying components. 



A. Example 1 

As the first example we consider a state defined as 

e W = C \0)d + C 1 \l) d + --- + C k - 2 \k-2) d + C k \k) d , (16) 

where Co,C k ^ and k is odd. Since k is odd ko = 1 
and Eq. (TTJ| shows that x 00 = C k ~\/C k = 0. Because 
P(£ N ;xoo) = Co ^ 0, £ cannot be a separable state. In 
a three-qubit case we see that, for example, 



£ 3 = C*o|000) + Ci|100) + C 2 |010) + C 3 |110) 
+ C 4 |001) +C 5 |101) + C 7 |111) 



(17) 



where C0C7 ^ cannot be separable. 

In order to compare our approach with other separa- 
bility tests, we now check the separability of £ 3 using an 
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alternative method. There exist various (partial) multi- 
partite separability criteria for mixed states (see, for ex- 
ample, While these are useful when mixed states 
are studied, in the case of pure states the most conve- 
nient separability check is usually the standard method 
of calculating the reduced single-qubit density matrices 
of the Af-qubit state. This view is supported by the fact 
that alternative pure state separability tests require ex- 
amining the properties of matrices that are higher dimen- 
sional than the two-by-two dimensional reduced single- 
qubit density matrices [l3|, [l8[ or require the calculation 
of the expectation values of operators expressed as tensor 
products of the Pauli spin matrices [HI, [Hi • This results 
in a complex calculation if a state that contains many 
freely varying components, such as £ 3 , is studied. For 
these reasons we now examine the separability of £ 3 by 
using the method of partial traces. Here and in what fol- 
lows we denote the reduced single-qubit density matrix 
pertaining to qubit j by pj. Now the indexing of qubits 
runs from to N — 1. The vector £ 3 is separable if and 
only if any two of the three density matrices po,pi, and 
P2 describe pure states. The state pj is pure if and only 
if det(pj) — 0, so if det(pj) ^ for at least one j, then 
£ is entangled. As an example we determine det(po). 
A simple calculation shows that the single-qubit reduced 
density matrix of the first qubit is 



Po 



\C \ 2 
C *Ci - 



-\c 2 \ 2 - 

csc 3 



|C 4 | 2 c ci + c 2 ci + c A c* 

-C|C 5 |Ci| 2 + |C 3 | 2 + |C 5 | 4 + |C 7 | 2 



(18) 

Using the inequality Re(C) < |C|, where C is an arbi- 
trary complex number, it can be shown that the following 
inequality holds for the determinant of po 

det(pn) > |C 7 | 2 (|C | 2 + |C 2 | 2 + |C 4 | 2 ) + (\C C 3 \ - |dC 2 |) 2 
+ (\CoC B \ - |dC 4 |) 2 + (\C 2 C 5 \ - \C 3 Ci\) 2 . 

(19) 

This is bounded below by |CoCV| 2 > 0, confirming the 
aforementioned result concerning the separability of £ 3 . 
Therefore a necessary condition for the separability of £ 3 
can be straightforwardly obtained using partial traces. 
However, the polynomial method provides a simpler sep- 
arability test in the present example. Even more so if 
instead of £ 3 the separability of the A^-qubit state £ N is 
studied. 



In the case of a three-qubit system this result means that 

e 3 = C o |000) + Ci|100) + C 3 |110) + C* 4 |001) 

+C 5 |101) + C 6 |011>, C 6 ^0, (20) 

cannot be a product state if Cq,Ci, or C3 is nonzero. 
If C4 = we have Xio = x\\ = 0, which means that 
all equal to zero. Then £ 3 cannot be separable 

unless all d except Cq are zero. The reduced single-qubit 
density matrices po,pi, and p 2 can be straightforwardly 
calculated and are not presented here. The determinant 
of p 2 is 

det( P2 ) = (|Co| 2 + |d| 2 + |C 3 | 2 )(|C 4 | 2 + |C 5 | 2 + |C 6 | 2 ) 
-\C Ct + CxCt\ 2 (21) 
> |C* 6 | 2 (|Co| 2 + |Ci| 2 + |C 3 | 2 ) 

+ |c 3 | 2 (|c 4 | 2 + |c 5 | 2 ) + (IOAI - I^QI) 2 , 

(22) 

where we have obtained a lower bound for the determi- 
nant in the same fashion as in the previous example. We 
reproduce the earlier result that £ 3 is necessarily entan- 
gled if C\, C 2 or C3 is nonzero. In order to determine the 
separability conditions in the case C4 = one has to cal- 
culate det(po) and repeat the above calculation for this 
quantity. The result agrees with the one obtained using 
the polynomial approach, that is, if £ 3 is separable and 
C4 = 0, then only Cq can be nonzero. We see that also 
in this case the polynomial approach provides an easier 
way to check the separability than the method of partial 
traces. 



C. Example 3 



As the final example we study a state given by 



E 



\i) d + e ie l 4 ^ 



i#0 : 4,8,...,2 w -4 



2"-l 



E i»% 



i=0 



1) l 4 ^- ( 23 ) 

i=0 



Now C(2N_ 1 )_ 2 j/C 2 iv„ 1 = 1 for all j, so Eq. (fTTj) gives 



B. Example 2 

In the second example we choose £ N such that the 
degree of P(£ N ;x) is k = 2 N - 2. Then kj = 1 - S oj . 
We assume that C 2 jv-i_ 2 (= C/j_ 2 iv-i) = 0, from which 
it follows that £(jv-i) m = for m — 0, 1, . . . , 2 7V ~ 1 — 
1. According to Condition (lb) the lowest order of the 
polynomial has to be at least 2 Ar_1 for the state to be a 
product state. Thus, if Ci ^ for at least one i such that 
< i < 2 N ' 1 , i ^ 2 N ~ 1 - 2, then £ N must be entangled. 



x jm =e 1 ^^, m = 0,l,...,2 J '-l, (24) 

where j — 0,1,2, ... ,N — 1. Using the sum formula of ge- 
ometric series we find that the characteristic polynomial 
can be written as 



x - 1 



x 1 ' - 1 



> x) = + (e 1} ^-r- (25) 

It is easy to see that for j = 2,3, . . . , N — 1 

P(Z N ;x jm ) = 0, m = 0,1,..., 2^-1, (26) 
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while 



P^;x 00 )=P(^;x 10 ) = P(^;x 11 ) 



>JV-2 



!)■ 

(27) 

= 27rn for 
shows that 



The state !; N is separable if and only if 8 
some integer n. If £ N is separable Eq. (fT3"|) 

e Ar = ®f=o 1 (|o)i + |i) J -). > 

Now the A~ reduced single-qubit density matrices of 
£ N can be straightforwardly determined. Lengthy calcu- 
lation shows that det(p ) = det(pi) = 2 2JV-3 (1 — cos#) 
and det(,Oj) = when j = 2, 3, . . . , N — 1, confirming 
the earlier result. In the present example the polynomial 
method does not seem to provide as obvious calculational 
simplification as in the previous two examples. 



D. Generalization to h-level systems 

We now briefly discuss a generalization of the separa- 
bility test to a system consisting of N copies of an /i-level 
system. We write the basis of a single /i-level system as 
{|0)/i, \ l)h, ... ,\h— l)h} and choose the basis vectors for 
the iV-partite system as \i)d = \ioii ■ ■ ■ iN-i)h, where 
ih J and ij <E {0, 1, 2, . . . , h — 1}. An arbitrary 



T N-1 



pure state can be expressed as 



h N -i 



c *\*)< 



(28) 



i=0 



Let 0s = 00^1 ' ' ' 0jv-i be a separable state where 0£ = 

a>j 1 0} h + bj 1 1} h + Cj 1 2) /( H h f/j | /i — 1 ) h . A straightforward 

calculation shows that 



3=0 



(.29) 



In order to establish a separability test, one has to express 
the roots of this polynomial in terms of the coefficients 
Co, Ci, . . . , C^jv_i. This is possible but complicated if 
2 < h < 6. If /i > 6, the roots cannot be calculated an- 
alytically and therefore cannot be written using the co- 
efficients Ci . Thus the separability test can be extended 
to systems containing less than six levels, but it is more 
complicated to apply than in the two-level case. An ex- 
tension is not feasible if the number of levels is equal to 
or larger than six. 



IV. NUMBER OF UNENTANGLED QUBITS 

Entangled states can be classified based on the number 
of unentangled one-qubit states. The state cj) is said to 
contain n unentangled qubits if it can be written as a 
product of n single-qubit states </>/ and an (N — n) -qubit 
state (j) n . In order to study the number of unentangled 
particles, we determine the characteristic polynomial of 
a state which separates as a product of a one-qubit state 



and an (N — l)-qubit state. We write <j) = <l>j<l> N ~ l , where 
4>j = aj\0)j + bj\l)j is the state of the qubit j and i^ -1 
gives the state of the rest of the qubits. As before, the 
degree of the polynomial is denoted by k. Using Eq. 
(|6]) we see that the characteristic polynomial of the basis 
states reads 



P(\i)d,; x) = P(\i ii ■ ■ ■ in-i)]x) 



x iq2 ^ixl x i 2 2 _ _ ^tiv-i^ 



(30) 



We write the (N — l)-qubit state as 
^n-i _y c ia ... ij _ l -i ]+1 ...i N _ 1 |*o * - * ij-iij+i ■ ■ ■ tN-i), 



so using Eq. (|30[) we find that 
P(<p;x) — bj(x — Xj m ) Ci ...i j _ 1 -i j+1 ...i N _ 1 



(31) 



X X 



i 2" + ---+i 3 _i2 J - 1 +i 3 + 1 2 J + 1 +---|-ijv-l2 J 



(32) 



where we have assumed that bj ^0, which is equivalent 
to kj = 1. We have also written (aj + bjX 23 ) = bj(— Xj m + 



-2^ 



Note that x 21 m is independent of to. If bj =0, we get 



an expression which is obtained by multiplying the sum of 
Eq. (|32|) by aj . Equation (|32|) shows that the polynomial 
P{4>; x)j[x v —xf m ) contains only those powers of x which 
do not have 2 3 in their binary representation and that 
Xj m is a root of P(<f>; x) for each to = 0, 1, . . . , 2 J — 1. 
Therefore, if kj = 1, necessary conditions for the qubit j 
to be unentangled with respect to the rest of the qubits 
are 

(Ila) P{4>] Xj m ) = for every to = 0, 1, . . . , 2° ' — 1. 
(lib) 2 J does not appear in the binary representations 
of the exponents of x in P((j>; x)/(x 21 — xf m ). 

If kj = there is only one condition, namely, 

(lie) 2 J does not appear in the binary representations 
of the exponents of x in P((f>; x). 

It is easy to see that these are also sufficient conditions. 
The number of unentangled qubits can be obtained by 
checking Conditions (Ila) and (lib) for every qubit j for 
which kj = 1 and Condition (lie) for the rest of the 
qubits. It is possible to extract information about the 
number of unentangled qubits without using Conditions 
(lib) and (He), namely, an upper bound for this quantity 
can be obtained by adding to the number of qubits for 
which (Ila) holds the number of indices j for which kj = 
0. This corresponds to assuming that either (lib) or (lie) 
holds for every qubit. 



A. Example 1 

As an example of the use of this method we consider 
the state given by Eq. (|2"3")l . Now k = 2 N — 1 and therefore 
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kj = 1 for every j. Equations ([26]) and (|27|) together with 
Condition (Ha) show that the number of unentangled 
qubits is at most N - 2 (N) if 9 ^ 2im (9 = 2vm). In 
order to simplify the polynomial P(£; x) we note that 

X ^ N - 1 = - l)( x 2 + l)( x 4 + l) . . . (x 2 ™" 1 + 1). (33) 

With the help of this and Eq. (|25l) we get 

P{£ N ; x) = (x + l)(x 2 + l)(x 4 + 1) • • • {x 2 ^ 1 + 1) 
+ (e w - l)(x 4 + l)(a; 8 + 1) • • • (x 2 "' 1 + 1) (34) 

Now x 2 ' — (xjm) 2 = x 2,3 + 1 when j > 2 and using the 
above equation one can see that Condition (lib) holds for 
j = 2,3,...,iV — 1 regardless of the value of 9. Further- 
more, (lib) holds for every j if 9 = 2irn. In conclusion, 
the qubits j = 2, 3, . . . , N — 1 are always unentangled 
with respect to the rest of the qubits and if 9 = 2m 
the state is separable. The same result can be obtained 
using the reduced single-qubit density matrices pj. The 
number of unentangled qubits is equal to the number of 
Pj for which det(pj) = 0. The values of these determi- 
nants have been presented in Example 3 and reproduce 
the aforementioned result. 

A necessary step in the calculation of the number of un- 
entangled qubits is to apply Condition (Ha) to all qubits 
j for which kj — 1. This is equivalent to checking the 
separability of the state. In addition to this, Conditions 
(lib) and (lie) have to be controlled. On the other hand, 
in the case of single-qubit reduced density matrices pj, 
the determination of the number of unentangled qubits 
does not require any additional operations in comparison 
with testing the separability. In both cases det(pj) has to 
be calculated. This suggests that the method of reduced 
single-qubit density matrices is preferable if the number 
of unentangled qubits is studied. 

V. CONCLUSIONS 

We have defined a mapping which associates pure N- 
qubit states with a polynomial. The roots of this polyno- 



mial determine the state completely and vice versa. The 
structure of the polynomial is inspired by the one used 
in the Majorana representation [lij [2(3] ■ The separability 
of a state can be studied by examining the properties of 
the roots of the corresponding polynomial. In particular, 
we have presented a method which establishes a neces- 
sary and sufficient condition for a given pure A^-qubit 
state cj) to be separable. This method provides a new 
point of view to the pure state separability and gives an 
alternative to the conventional separability test of cal- 
culating the reduced single-qubit density matrices of the 
state. The separability of </> can be determined by check- 
ing whether the numbers Xj m , defined in equation (|14[) . 
are roots of the polynomial P(4>', x) of equation ((2J. Both 
the numbers Xj m and the polynomial P(<j>; x) can be eas- 
ily obtained as a function of the components of the state 
<p. We have illustrated through examples that in some 
cases the polynomial separability test is easier and faster 
to apply than the method of reduced single-qubit den- 
sity matrices. We have also shown how the number of 
unentangled qubits can be obtained with the help of the 
polynomial P{<p; x). It seems, however, that for this task 
the method of single-qubit density matrices is preferable. 
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